rvis 13

Load Packages and Data

Daten angucken

str(death)

'data.frame':   250 obs. of  4 variables:
 $ count   : int  3 2 1 1 4 2 2 2 3 2 ...
 $ geometry: chr  "<Point><coordinates>-0.13793,51.513418</coordinates></Point>" "<Point><coordinates>-0.137883,51.513361</coordinates></Point>" "<Point><coordinates>-0.137853,51.513317</coordinates></Point>" "<Point><coordinates>-0.137812,51.513262</coordinates></Point>" ...
 $ lon     : num  -0.138 -0.138 -0.138 -0.138 -0.138 ...
 $ lat     : num  51.5 51.5 51.5 51.5 51.5 ...

str(pump)

'data.frame':   8 obs. of  4 variables:
 $ count   : int  -999 -999 -999 -999 -999 -999 -999 -999
 $ geometry: chr  "<Point><coordinates>-0.136668,51.513341</coordinates></Point>" "<Point><coordinates>-0.139586,51.513876</coordinates></Point>" "<Point><coordinates>-0.139671,51.514906</coordinates></Point>" "<Point><coordinates>-0.13163,51.512354</coordinates></Point>" ...
 $ lon     : num  -0.137 -0.14 -0.14 -0.132 -0.134 ...
 $ lat     : num  51.5 51.5 51.5 51.5 51.5 ...

Get Maps

berlin.map <- get_map("Mitte, Berlin",
                        zoom = 16, 
                        maptype = "terrain", 
                        source = "stamen")

anzeigen lassen

# mit cowplot
plot_grid(plotlist = plots, ncol = 2, labels = c("Google Maps", "Stamen Maps", "Berlin Map"))

Fehler in plot_grid(plotlist = plots, ncol = 2, labels = c("Google Maps",  : 
  Objekt 'plots' nicht gefunden

weitere Layer hinzufügen

Adressen Geokodieren

home <- geocode(location = "Birkenstr. 11, Berlin", source = "google")

Source : https://maps.googleapis.com/maps/api/geocode/json?address=Birkenstr.%2011%2C%20Berlin

Adresse zu Map hinzufügen

Shapefiles

bb.shp <- readOGR(dsn = "data/envvk9q3g/DE_BB_Mroad_Lnight.shp", layer = "DE_BB_Mroad_Lnight")

OGR data source with driver: ESRI Shapefile 
Source: "data/envvk9q3g/DE_BB_Mroad_Lnight.shp", layer: "DE_BB_Mroad_Lnight"
with 8117 features
It has 6 fields

1.Feb Maps II

germany <- get_map(location = "Germany",
                   zoom = 6,
                   maptype = "toner",
                   source = "stamen")

Source : https://maps.googleapis.com/maps/api/staticmap?center=Germany&zoom=6&size=640x640&scale=2&maptype=terrain
Source : https://maps.googleapis.com/maps/api/geocode/json?address=Germany
geocode failed with status OVER_QUERY_LIMIT, location = "Germany"Fehler in data.frame(ll.lat = ll[1], ll.lon = ll[2], ur.lat = ur[1], ur.lon = ur[2]) : 
  arguments imply differing number of rows: 0, 1

Have a first look

Katre laden und als kontext Layer verwenden

# Stadtnamen ergänzen
city[c(3, 5, 6), "Hauptstadt"] <- c("Berlin", "Bremen", "Hamburg")

Fehler in city[c(3, 5, 6), "Hauptstadt"] <- c("Berlin", "Bremen", "Hamburg") : 
  Objekt 'city' nicht gefunden

Daten räumlich machen

# General logic for geocoding
geocode(location = "Stuttgart", output = "latlon", source = "google", inject = ", Germany")

Source : https://maps.googleapis.com/maps/api/geocode/json?address=Stuttgart&,%20Germany

Für viele Städte

# for loop with if-condition for warnings
# run several times until no more NA's are returned
for(i in 1:nrow(city.df)) {
  if (is.na(city.df$lon[i])) {
    result <- tryCatch(geocode(city.df$Hauptstadt[i], output = "latlon", source = "google", inject = ", Germany"),
                       warning = function(x) data.frame(lon = NA, lat = NA))
    city.df$lon[i] <- as.numeric(result[1])
    city.df$lat[i] <- as.numeric(result[2])
  }
  i <- i + 1
  Sys.sleep(sample(seq(.5, 2, 0.5), 1))  
}

Source : https://maps.googleapis.com/maps/api/geocode/json?address=Stuttgart&,%20Germany
Source : https://maps.googleapis.com/maps/api/geocode/json?address=M%C3%BCnchen&, Germany
Source : https://maps.googleapis.com/maps/api/geocode/json?address=Berlin&,%20Germany
Source : https://maps.googleapis.com/maps/api/geocode/json?address=Potsdam&,%20Germany
Source : https://maps.googleapis.com/maps/api/geocode/json?address=Bremen&,%20Germany
Source : https://maps.googleapis.com/maps/api/geocode/json?address=Hamburg&,%20Germany
Source : https://maps.googleapis.com/maps/api/geocode/json?address=Wiesbaden&,%20Germany
Source : https://maps.googleapis.com/maps/api/geocode/json?address=Schwerin&,%20Germany
Source : https://maps.googleapis.com/maps/api/geocode/json?address=Hannover&,%20Germany
Source : https://maps.googleapis.com/maps/api/geocode/json?address=D%C3%BCsseldorf&, Germany
Source : https://maps.googleapis.com/maps/api/geocode/json?address=Mainz&,%20Germany
Source : https://maps.googleapis.com/maps/api/geocode/json?address=Saarbr%C3%BCcken&, Germany
Source : https://maps.googleapis.com/maps/api/geocode/json?address=Dresden&,%20Germany
Source : https://maps.googleapis.com/maps/api/geocode/json?address=Magdeburg&,%20Germany
Source : https://maps.googleapis.com/maps/api/geocode/json?address=Kiel&,%20Germany
Source : https://maps.googleapis.com/maps/api/geocode/json?address=Erfurt&,%20Germany

Ausläneranteil darstellen

ger.map +
  geom_point(data = city.df, mapping = aes(x = lon, y = lat, size = Ausländeranteil), color = "gold2", alpha = .5)

Fehler in eval(expr, envir, enclos) : 
  Objekt 'Ausländeranteil' nicht gefunden

geht auch mit ggmap

bl.shp <- readOGR(dsn = "data/VG250_BL_2017",
                  use_iconv = TRUE, encoding = "UTF-8", stringsAsFactors = FALSE)

OGR data source with driver: ESRI Shapefile 
Source: "data/VG250_BL_2017", layer: "VG250_Bundesländergrenzen_2017"
with 35 features
It has 26 fields
Integer64 fields read as strings:  OBJECTID ADE GF BSG IBZ 

cleaning

# (2) ggfortify
bl.for <- fortify(bl.shp)

Regions defined for each Polygons

plot it

Shapefile angucken

using tmap

tm_shape(shp = bl.shp) +
  tm_polygons()

mergen

bl.shp <- sp::merge(bl.shp, city.df, by.x = "GEN", by.y = "Land")

Fehler in sp::merge(bl.shp, city.df, by.x = "GEN", by.y = "Land") : 
  Objekt 'city.df' nicht gefunden

TMap stellt daten zur Verfügung WHOOOOOP